home *** CD-ROM | disk | FTP | other *** search
/ Mac Mania 5 / MacMania 5.toast / / Tools&Utilities / Plotfoil 3.2 / NACA-info < prev    next >
Encoding:
Text File  |  1995-09-09  |  2.3 KB  |  59 lines  |  [TEXT/ttxt]
  1.  
  2. /*
  3.  * naca.c
  4.  *
  5.  * Generate points for the NACA 4 and 5 digit airfoils.
  6.  *
  7.  *                        NACA 4-digit Airfoils
  8.  *
  9.  * The first digit gives the camber (in % chord); the second the x-value
  10.  * of the max. camber point (in tenths chord); the last two are the thickness
  11.  * (in % chord.). Thus a NACA 4312 has a 4% camber, with max. camber at 30%
  12.  * chord, and has a thickness of 12%.
  13.  *
  14.  * The camber line consists of two parabolae that are joined smoothly with
  15.  * horizontal slope at the max. camber point (x_a, y_a). Thus the equations
  16.  * for the camber line become:
  17.  *
  18.  *    y = y_a/x_a^2 x (2x_a - x)                   0 <= x <= x_a
  19.  *    y = y_a / (1 - x_a)^2 (1 - x)(1 + x - 2x_a)  x_a <= x <= 1
  20.  *
  21.  * This camber line deforms the standard NACA streamline shape:
  22.  *    +-d = t [ 1.4845 sqrt(x) - 0.63 x - 1.758 x^2 + 1.4215 x^3 - 0.5075 x^4 ]
  23.  *
  24.  * So for each point x, find the point on the camber line c(x) and the normal
  25.  * to it. Then along the normal find points at a distance of +- d, the
  26.  * thickness function. (This is an approximation, for small camber.)
  27.  *
  28.  *
  29.  *                NACA 5-digit Airfoils
  30.  *
  31.  * The 4-digit airfoils have the property that at the max. camber point the
  32.  * radius of curvature changes discontinuously. This implies a discontinuous
  33.  * centripetal acceleration in the airstream; to ameliorate this, the camber
  34.  * line is modified to one of two forms:
  35.  *
  36.  * i)  the "simple" camber line: this consists of a cubic joined to a straight
  37.  *     line. At the joining point, both the slope and curvature are continuous.
  38.  * ii) the "reflexed" camber line: this consists of two cubics joined together
  39.  *
  40.  * The second digit is the design lift coefficient.
  41.  * The third digit is either 0 (simple camber) or 1 (reflexed camber). The
  42.  * last two digits are thickness just as in the 4-digit case. The first two
  43.  * digits define the camber line. The second digit is multiplied by 5 to give
  44.  * the position of max camber. The max. camber is defined according to this
  45.  * table:
  46.  *  second digit ->     1    2    3    4    5
  47.  *  first digit
  48.  *        2    1.1    1.5    1.8    2.1    2.3
  49.  *        3        2.3    2.8    3.1
  50.  *        4        3.1    3.7    4.2
  51.  *        6        4.6    5.5    6.2
  52.  *
  53.  * Unfortunately I don't know how to generate the camber line from this
  54.  * information; for each cubic now we have two points and one slope, which is
  55.  * not enough to uniqely determine the curves.
  56.  *
  57.  */
  58.  
  59.